Cas 
hi hdits tes 
es 


Gx wpris 
UNIOMASTTATIS 
AABERTAEDSIS 


SPECIAL COLLECTIONS 


UNIVERSITY OF ALBERTA LIBRARY 


REQUEST FOR DUPLICATION 


I wish a photocopy of the thesis by 


» 2, 
\ on {2 ot 


(author ) 


5 A, 
1 


Aactics in th hortyt Arch 


entitled. Chat 


The copy is for the sole purpose of private scholarly or scientific study 
and research. I will not reproduce, sell or distribute the copy I request, 
and I will not copy any substantial part of it in my own work without per- 
mission of the copyright owner. I understand that the Library performs 
the service of copying at my request, and I assume all copyright responsi- 
bility for the item requested. 


Digitized by the Internet Archive 
in 2023 with funding from 
University of Alberta Library 


https://archive.org/details/Pelot1982 


- Ly rea ( i 


re > ‘ Ve y @ “7; 
om i a as = - 
. r 7 ee 


BL 
a 


; ; 
4 


7 


Eo EUN DVR Row Ye OP AL ar Ro A 


RELEASE FORM 


NAME OF AUTHOR Ron Pelot 
TITEBSOr THESIS Flow Characteristics in the Aortic Arch 
DEGREE FOR WHICH THESIS WAS PRESENTED Master of Science 
YEAR THIS DEGREE GRANTED SVolwa ieley, fekst2 
| Permission is hereby granted to THE UNIVERSITY OF 
ALBERTA LIBRARY to reproduce single copies of this the- 
gis and to lend or sell* such copies for private, scho- 
larly or scientific research purposes only. 
The author reserves other publication rights, and 
neither the thesis nor extensive extracts from it may 
be printed or otherwise reproduced without the author's 


written permission. 


© SrasTrwes Ww" om tee 4 ese: = halael® a 
=r 8. er “a @: — oFoe¥e eee ne ee) ae 


a 


& 
4 


THE UNIVERSITY OF ALBERTA 


Flow Characteristics in the Aortic Arch 
by 
ee 
fof WP eRonePelot 


4 ; f 
~. J 
Na ae 


A THESIS 
SUBMITTED TO THE FACULTY OF GRADUATE STUDIES AND RESEARCH 
IN PARTIAL FULFILMENT OF THE REQUIREMENTS FOR THE DEGREE 


OF Master of Science 


Mechanical Engineering 


EDMONTON, ALBERTA 


Spring, 1982 


mee 


Ag 


etéanT. 4 j 
spade 2 stat eh 7 anieniaed 


THIMITMON IH fa 10 TAP patle jects V2 


| 


a°; ot Gw a bs ase 2 


THE UNIVERSITY OF ALBERTA 


FACULTY OF GRADUATE STUDIES AND RESEARCH 


The undersigned certify that they have read, and 
recommend to the Faculty of Graduate Studies and Research, 
for acceptance, a thesis entitled Flow Characteristics in 
the Aortic Arch submitted by Ron Pelot in partial fulfilment 


of the requirements for the degree of Master of Science. 


Jvecep of to ealp. Se 


7 7 a - —— rn 


7 
- 


_ -) 7 
+) a 
7 _ os 


Abstract 


The flow patterns in a model of the aortic arch are studied 
by visualization under steady and pulsatile conditions. Par- 
ticular attention is paid to secondary flow patterns, 
separation zones within the arch and flow distribution as a 
function of Reynolds number. 

The flow is characterized by secondary flow vortices 
induced by centrifugal forces, a core of fluid which remains 
potential through the length of the arch and backflow areas 
along the inner curve wall and outer curve wall. The inner 
curve separation zone can be detected when the Reynolds num- 
ber exceeds one thousand during steady flow, but persists 
throughout the period for pulsatile tests. 

The sites of backflow coincide with the localization of 
atherogenesis and provide further evidence in establishing a 
correlation between hemodynamics and disease development in 


the aorta. 
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favintroduction 
Much s7o fo@rtne progress in medicine has come from exchanging 
acute medical problems for chronic ones. The increasing abi- 
lity to understand, diagnose and cure diseases and defects 
has relegated many previously '‘incurable' afflictions to the 
realm of improbability. The human cardiovascular system 
however remains a major cause of disability and death. 

Engineering research in the field of hemodynamics (hy- 
drodynamics of blood flow) provides information about the 
characteristics of blood flow through arteries and veins. 
Coupled with medical research, hemodynamics leads to a_ bet- 
ter understanding and possible correction of physiological 
impairments. 

One such problem concerns the development of athero- 
sclerosis in the arterial system. These lesions inhibit the 
flow and can ultimately block an artery to the point where 
the blood flowrate 15 critically low. There are several 
theories which have been proposed to explain lesion forma- 
tion. Most of them have a biochemical basis but none com- 
prise a complete explanation of this complex phenomenon. 

The atherosclerotic lesions usually originate at spe- 
cific sites in the large arteries, particularly at bends and 
branches. Therefore, hemodynamic theories have been advanced 
pPOerelatewlesion formationsto the dynamic sbehaviour of blood 
abethese Sites. The aorta, an inverted U-shaped artery ori- 
Ginatincmec CC ICM Le hLVentiiCLem(N1GQ.iW my m@elomOlLeDaleicular 


interest since it encompasses several potential lesion 
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Sites. 
A complete picture of the flow patterns in large 
arteries could provide clues to the cause and formation of 


atherosclerotic lesions. 
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Figure 1.1 Human heart and aorta 
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2. History of the problem 

MNeseXac te nauIrenotetnesvel Ocuty: d1Stripumvonee nes Lhesmaonca 
plays an important role in the study of atherosclerotic for- 
mations. There are two major opposing theories which attempt 
to correlate primary atherosclerotic lesions with the 
dynamic behaviour of blood. Fry(1968) demonstrates a large 
increase in endothelial damage at high shear rates and sub- 
sequent deposition of blood elements at these locations. 
Conversely, Caron fete ale Ho Oem conmencswamLnatemmncip vent 
atheroma (atherogenesis) coincides with low sanguinary shear 
Gale Sam which fans to Seema ent scour the arterial walls and 
prevent excessive plaque build-up. Scarton et al.(1976) 
postulates that one theory does not necessarily preclude the 
other, and that either effect may be the causal agent of 
atherogenesis at various sites. 

The distribution of atherosclerotic lesions during the 
various stages of their development has been documented by 
Several researchers. Some studies have been conducted on 
animals (ex. Rodkiewicz,1975), others on human cadavers 
(Shah et al.,1976) and some information apse severe condi- 
tions has been gleaned by electrocardiographic and other 
diagnostic procedures (Rushmer,1970). Gresham( 1976) presents 
an excellent review of experimentally induced and spontane- 
ous atherosclerosis, particularly in primates. 

Although many of the observations conflict, the 
affected areas in the aortic arch region generally include 


the inner curve of the arch and less severe formations on 
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ENCROUCCKEDORtwiOn Oia theacurve awhile them Stcew walls. are 
relatively unaffected. As well, there appears to be some 
build-up at the entrances of the large branches and on _ the 
proximal walls of these branches just beyond their inlet. 

In light of these obServations, a fluid flow investiga- 
tion has taken place on several fronts. One approach con- 
Sists of an analytical and/or numerical analysis making use 
of the Navier-Stokes equations with several simplifying 
assumptions. In its most basic form, the aortic arch can be 
modeled as a curved rigid tube. Dean(1927) considered this 
problem for steady laminar flow and demonstrated the nature 
of the secondary flow induced by centrifugal effects. The 
nature of pulsatile flow in circular straight pipes was exa- 
mined by Womersley(1955a) and Uchida(1956). Unfortunately, 
the integration of these two analyses (Dean's and 
Womersley's) has defied rigorous and complete treatment due 
to its mathematical complexities. Smith(1975) has considered 
pulsatile flow in curved tubes, Chandran et al.(1974) ana- 
lysed flow in curved elastic tubes and Singh(1974) looked at 
entry flow in a curved pipe. All of these mathematical 
treatments rely however on numerous simplifying assumptions 
mcludind sconstante diameter pe slardemscurvaturemErabiogsand 
absence of branches. Curvature ratio iS a non-dimensional 
parameter defined as the radius of curvature of a curved 
tube divided by the internal radius of the tube. 

A second line of research involves direct measurement 
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blood flow in vivo have proven to be very difficult. The 
Slow development of non-invasive probes had precluded most 
investigations on healthy humans) where the term healthy 
here denotes any person without significant heart diseases 
or defects. The time-varying contours of intravascular and 
intracardiac pressures have been easily measured externally 
and extensive investigations of these quantities have been 
accomplished by non-invasive techniques. By approximating 
anterval’itoweas pulsatile= «flow through@yinfimite, ainigid® 
circular tubes, Womersley(1955b) mathematically derived the 
velocity curve from the pressure measurement. However, the 
many Simplifications inherent in applying such an analysis 
torthe aorta prompted “further developments in meaSuring the 
time-varying velocity more directly. 

Several instruments were developed and improved over 
the years, initially focussing primarily on periarterial 
coupling. The electromagnetic flowmeter (Kolin, 1936,1937), 
the pulsed ultrasonic flowmeter (Franklin,et al.,1959), the 
Doppler flowmeter (Franklin et al.,1961) and the hot-wire 
anemometer (Machella,1936; Ling et al.,1968) have all been 
considerably refined since their conception. These flowme- 
ters have a common major drawback in that they require expo- 
sure of the artery. Therefore, their application has been 
limited to animals and to humans undergoing surgery. 
Although they provide an abundance of valuable information, 
these instruments again do not produce the velocity trace of 


a healthy human. 


7iJo Li Gai |! 


ty aid Pon va 


a 


Miniature versions of these flowmeters and advances in 
Surgical techniques have permitted the introduction of such 
devices into arteries by means of a catheter (Mills,1966; 
Metis andeShillingftord); 1967). The instantaneous velocity at 
various sites in healthy humans has thus been recorded by 
several researchers using intravascular probes (Gabe 
Sema ln 1969 MeMitllsmeetmecte w1970- =Nichololetura lemoy7a a 
However, these instruments initially could not measure spa- 
Dia kevanvatromeolevelocityreacross thertubes  shneim applica— 
tion assumes uniform velocity across the lumen 
(McDonald,1974) and gives no indication of secondary flow 
patterns. This information is required for a complete map- 
ping of the flow field. 

| Fortunately, further advances in instrumentation has 
Stimulated more studies which yielded some velocity profiles 
Pomanimameareeriess  Seedmand wood ( 1970, 107 |) stransversed: the 
aortae of dogs with a hot-wire anemometer by directly punc- 
turing the vessel wall. Farthing and Peronneau(1979) applied 
pulsed-Doppler ultrasonic transducers to the exterior of 
Canine aortae which yielded velocity profiles at various 
Sites and at different times throughout the cycle. Although 
the morphologies of the canine and human aortae are somewhat 
different, these results provide invaluable insights into 
the flow behaviour. Neither of these techniques can detect 

Separation zones nor discern secondary flow patterns. 
A final way to study the flow patterns involves the use 


of physical models. A drawback of this approach concerns the 
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limited degree of simulation obtainable in the laboratory, 
whether it be morphological reproduction, dynamic similarity 
Or correspondence of various other characteristics. The 
great” advantage’ of this method lies in’ the ability? to con- 
trol parameters at will and study their individual effects. 
A number of researchers have considered the aortic arch 
flow from several points of view. For example, Rodkiewicz 
et al.(1976a) examined separation zones in a plastic model 
of the human aorta under steady and pulsatile flow condi- 
tions using hydrogen-bubble visualization. Austin and Sead- 
er(1973) experimentally studied the entry region of a curved 
tube while Scarton et al.(1976) revealed flow details ina 
curved glass tube. Since published results do not always 
agree due to differing experimental conditions, parameters 
and assumptions, this line of research’ =continues  untrie a 


Single clear pattern emerges. 
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3. Background 
A complete description of the aortic flow has not yet been 
formulated but several individual aspects have been revealed 
in various studies. Those characteristics which have been 
reasonably well established are reviewed in this section. 
Throughout this document, the terms mean and average are not 


used synonymously. They are respectively defined as: 
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faverage [t=cst Sy] 5 


where £ is anrarbitrary function of space(s) and time(t). 

Mode Wii numthemaont (cmb lLoodmrlowlwicheamsteady emt lowrate 
can be justified by several considerations. The pulsatile 
flow component is of the same order of magnitude as the mean 
forward flow. However, spatial variations in the wall shear 
stress depend primarily on the boundary layer growth due to 
the steady flow component. The fluctuating shear stresses 
Gausedsby sthesoscillatony component ane much jlessmestnongly 
dependant on position. Therefore, the mean shear stress has 
a spatial distribution similar to that of steady flow (Caro 
et al.,1971). 

Likewise, Rodkiewicz et al.(1976a) has demonstrated 
that the locations of the separation regions can be similar 


for steady and pulsatile flow, notwithstanding their size 
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variation under hluceuating conditionse mTheresimilarity 
between the flattened velocity profile in the core of oscil- 
latory flow (Atabek et al.,1964) and the flattened inlet 
velocity profile of steady flow further admits constant flow 
Simulation. For these reasons, a steady flow model can 
demonstrate many characteristics of pulsatile flow. 

Since the left ventricle acts as a reservoir and the 
aortic valve opens and closes very quickly, the velocity 
profile at the inlet to the ascending aorta is reasonably 
flat. This uniformity is further insured by the presence of 
the sinus of Vasalva, a deep narrow expansion of the aorta 
at the level of the coronary arteries which tends to absorb 
any jet effects during systole and also damps some vorticity 
Obiginatingpatwenes valvesn(McDonald 4974) ee (rigsc.a). —rhus 
description has been supported by measurements in dogs by 
Seed and Wood(1971) and Farthing and Peronneau(1979) which 
Show a remarkably flat profile in the ascending aorta. 

As mentioned previously, flow in the aortic arch is 
Similar in some respects to flow through a curved tube. In 
potential flow through a curved tube with a uniform velocity 
profile at the entrance, the total head pt+tpv* remains con- 
Stant along any streamline assuming the centreline lies in 
the horizontal plane. The static pressure increases with 
radius. from the centre of curvature to balance the centri- 
fugal force as the flow passes into the bend. Consequently 
the velocity increases at the inside of the bend while at 


the outer curve some velocity head tpv? is converted to 
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pressure head 'p' (Miller,1971). Coming out of the bend, the 
COnVeusSemocCursecCO recat imosunlt OLMMmpLornvemml gy Gk 20s 

For fully developed flow of a real fluid entering a 
Pipe bend, the slower moving fluid near the walls is acted 
Onmeipyana “Lowerecentriftugalmtorce than iiastem movingellurdy in 
the core. This central fluid is displaced outwards and the 
Slow-moving boundary fluid moves inwards around the walls 
towards the inside of the bend describing a double helix 
(Fig.3.3). Dean(1927) theoretically analysed this problem 
for laminar flow while White(1929) experimentally investi- 
gated the flow behaviour as reviewed by Schlichting(1968). 

Although the flow is much more complex for ae real 
fluid, pressure differentials between the inside and outside 
walls differ little from those calculated assuming an ideal 
SVU SOA tMERENe® OCUbLEt MES rOMEaebend nat nem: LOW mOnmte he =ans ice 
experiences a pressure gradient of the same magnitude as 
Ehate "predictedy ior. angrdeatetlurdieligsumiicientmveloci ty 
head is available to be converted to pressure head, the flow 
can negotiate the rise in pressure, otherwise it separates 
from the inner wall. 

Since the fluid entering the aorta has a reasonably 
“uniform vellocityapnotrie, ahhesintluencerotgcurvaturesiseless 
than in fully developed flow. This flow therefore does not 
fall into either of the above categories but rather exhibits 
Gharacteristics of booths Theoretically va potential core 
Glowesnouldmpersistetheougner nem length ee Oluee hese arch sine 


entry arc required for the flow to fully develop in a curved 
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bupewase quoted by Scartonmetal., (197/6)mismoiven by: 


a) Re (2) 


where @ is degrees of arc meaSured from the entrance. For 
Re=1000 and (R/c)=3/4, fully developed flow would not be 
established until 9=424°. Since the aorta describes an arc 
of only 180° approximately, the flow is always developing. 
Since the core remains potential, one would expect that 
the velocity at a cross-section varies inversely with the 
radius of curvature, which is supported by Farthing and 
Peronneau's(1979) measurements in dogs. However, the gradu- 
ally developing flow will assume some characteristics of 
real fluid in curves and helices can thus form along the 
boundaries (Rodkiewicz et al.,1976a; Timm,1942). 
Rodkiewicz(1976b) also gives evidence of separation and 
Stagnation regions. Shear stresses act away RRonmecum StaGias 
EL Onewarces (Pique. 4ajewhiichecan defornmthemwallthuse initia 
ting a separation zone (Fig.3.4b). Conversely, shear 
stresses compress the wall lining towards a separation arc 
(Fig.3.5a) and as the wall loses resiliency a bulge may form 
thereby enlarging the separation region (Fig.3.5b). These 
two effects could intensify atherogenesis by particle circu- 
lation and deposition rather than removal and by damaging 
the protective biochemical defences of the wall cells 
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Figure 3.1 Sinus of Vasalva 
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Figure 3.3 Fully developed flow(after McDonald, 1974) 
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Figure 3.4 Stresses induced by stagnation(Rodkiewicz, 1980) 
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Figure 3.5 Stresses induced by separation(Rodkiewicz, 1980) 


4. Experimental Parameters 
From one individual to the next, there is evidently a great 
variation of the parameters which are related to flow in the 
aorta. However, a good starting point when dealing with 
experimental models is to set each parameter at the average 
value found in nature, and then in subsequent studies exa- 
mine the effects of varying each parameter within its full 
physiological range. Unfortunately, large discrepancies in 
these parametric values are found in literature, therefore a 
comprehensive review is presented here. A summary of the 
average values used in this experiment as well as the normal 


physiological range is given in Table 4.1. 


4.1 Geometric Parameters 


4.1.1 Morphology 

The aorta arises from the left ventricle of the heart, 
passes upwards (ascending aorta), bends over (aortic arch) 
then passes down through the thorax and abdomen (descending 
GOREa) PCG rants) 195) Grin lr ne Me InNOG elma TemE tne Sans EUG Y 
represents the aorta from the root to the proximal part of 
the descending aorta. The arch is characterized by continu- 
ous convergence and by three large branches originating at 
the top of the curve. There are two carotid arteries which 
Supplysmostect the blocd tothe headtandwbrain = Theslert and 


right subclavian each lead to an arm (or foreleg in the case 
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Parameter 


Radius 


Density 


Dan MebGtehe 


Kite SC. 


Flowrate 
Velocity 


Reynolds 
number 


PeeeOQuUenc y 
parameter 


Table 4.1 Experimental Parameters 


Symbol 


Re 


Average 
value 


1.45cm. 
feOSeovcmn: 
de By SIONDE 

OY SONG SNE 
6eltoy/ min 
15.1em/sec 


1000 


18.8 


Normal 
range 


ORGasleocm. 
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2h Bie 1 SIO SRE + 

4 ees easy Seabiey Ane 
5.1-25cm/sec 


SO 0s 000 


10-30 


ifs) 


of animals). The first branch from the aorta is the brachio- 
cephalic artery which bifurcates into the right carotid and 
right subclavian arteries. For convenience, these branches 
and the descending aorta have been labelled 1A,1,2,3,4 and 5 
respectively (Fig.4.2). The arch centerline does not lie 
completely in one plane but 'twists' backwards by approxim- 
ately a 15° angle (Fig.4.3). 

For this project, the coronary arteries are " omitted 
Since they are difficult to replicate faithfully and do not 
Significantly affect the flow patterns in the arch. The 
ratio of their combined areas to that of the aortic inlet is 
less than 4% and together they consume only 3-4% of the 
blood flowing from the left ventricle. They originate at the 
Sinus of Vasalva which almost completely damps any disturb- 
ing effects they might have on the downstream flow patterns. 
The twist is also neglected, hence the model centreline lies 
in one plane. The effect of twist on the flow pattern is 


left as another study. 


4.1.2 Aortic Dimensions 

The inside diameter at the base of the human aorta is 
largely underestimated in literature. An average value of 
2.0-2.5cm. is often quoted and although this certainly falls 
within normal physiological range it is not representative 
of the mean. This discrepancy could arise from taking mea- 
Surements from cadavers where the tissue has inevitably 


contracted or by estimations from X-ray photographs where 
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the boundaries are quite indistinct and magnification errors 
Ma yeOCCurs 

Dotter and Steinberg(1949) angiocardiographically mea- 
Sured one hundred normal aortae and obtained a mean of 
28.6mm. with a range of 16-38mm. for the inside diameter at 
the mid-ascending aorta. Furukawa et al.(1976) found similar 
values from echocardiograms of 194 normal aortae. Their mea- 
Surements yielded a mean internal diameter of the aortic 
POCtmeat eC Hem enGunOlmmalaStoLle» -OlT 2 Ommem—mewitneasrange so 
14-39mm. 

Since the present model originates between the aortic 
root and the mid-ascending aorta, and since the diameter is 
2-3% larger during systole, a value of 29mm. has been 
adopted as a suitable average. 

The dimensions of the aortic arch and its derivative 
branches are’ shown in Figure 4.4. The data have been culled 
from several sources (Reich,1949; Meschan,1975; Reul 
Creal 1Ogtew Rodkiewiczeetwal 10/79) dues cvomaslack ol mmaqree— 
ment on average values and absence of a single comprehensive 


review. 


4,2 Flow Distribution 

The volumetric distribution of the flow in each branch 
Tomanweimpontant. factor anauenesreplicatron Omethe=actualst low 
Parcenns mmeACCOrding, = cOsmReulmer ae S7ayeedboutet%eor the 
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subclavian artery, while approximately 7% is diverted 
CNroOugnmwaeachmecarotldmmartery sleéavingamysoseco continue 1s 
course through the descending aorta. This agrees with values 
Gited »yin. -Guyton(1956)" sand, Ruch and. Patton(1973) although 
these sources do not provide a complete list. The model in 
this experiment has therefore been calibrated to the appro- 


priate distribution(Fig.4.2). 


4.3 Fluid Properties 

The dynamic viscosity of blood has been one of the most 
critically and comprehensively studied blood flow parameters 
(ex.Chien,1979). Blood is essentially a viscoelastic fluid 
and at a given temperature the fundamental factors affecting 
its viscosity are cell concentration, plasma viscosity, cell 
deformation and cell aggregation. However, for a given hema- 
tocrit (% cells in plasma), at shear rates between 107-10? 
sec”™', in tubes with internal diameter greater than I1mn., 
the viscosity becomes virtually constant (McDonald, 1974; 
Haynes,1961; Bayliss,1962). The lower limit on tube size is 
required since tests have shown that as the ratio of blood 
cell diameter to tube diameter becomes small the viscosity 
GendSeetowards, sa ~CONnStant value, all other things being 
equal. Since the aorta satisfies the stipulated criteria and 
hematocrit varies little between healthy individuals, it is 
possible to assume a constant viscosity thereby allowing the 


use of any Newtonian fluid (such as water) as the working 
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medium. 

Rosenblatt(1965) quotes a mean viscosity jy =4.59cP. at 
body temperature (37°C) with a standard deviation of o =.485 
for 117 normal male adults studied. The average value for 
females is about 14% lower but since the occurence of severe 
atherosclerosis predominantly affects men, this study will 
adopt values representative of adult males, ergo U =4.6cP. 

In ponpaen conten iod density has a much smaller physio- 
logical range and is much more amenable to accurate measure- 
ment. Most sources (Muirhead et al. ,1946; Gray and 
OLAn@he, WSvise))) Concun on an average blood density of 
Das Oo0g7 cmaawi thima ieee range 1.052<p<1.064 g/cm?. 

A proper evaluation of these fluid properties is crit- 
ical to ensure dynamic similarity of the actual aorta and 


the model (Section 4.8). 


4,4 Flowrate 

The variability of the average flowrate in normal phy- 
Siological situations exceeds that of all the other relevant 
parameters. Eclipsing the normal variation from one indivi- 
dual to another, the large range is primarily a function of 
activity level as well as sex, physical size, age and other 
“minor factors (Detweiler,1972). Although the range extends 
Promises Dreres/min. 8£0m a ssmallw adult woman eeaty rest esto 
45 litres/min. for a well-conditioned athlete engaged in 


vigorous exercise, one may assume that only a small range is 
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applicable to the study of atherosclerosis. Since 
middle-aged men and elderly people are most susceptible to 
the debilitating effects induced by advanced atheroma, and 
many have a low level of physical activity, a value of 
6 litres/min. is deemed representative of resting, sitting 
and light walking conditions (Guyton,1956; Bergel,1972) and 
a range from 4-8 litres/min. would encompass a large percen- 
tage of normal subjects. The accuracy of this assumption is 
not critical however since the overall flow field should not 


be affected significantly by this parameter. 


4.5 Velocity 

Values of velocity in the aorta found in literature are 
difficult to compare since they are alternately given as 
peak value, mean systolic or mean velocity. They are also 
estimated at various cross-sections in the aorta and a uni- 
form velocity profile is usually assumed at the site in 
question. The most consistent and expedient approach is to 
calculate the mean velocity from the flowrate and ascending 


aorta dimensions: 


_ 6.0 L/min. 15.1 cm/sec (4.1) 
alt) S)Iran 


This estimate falls somewhat below commonly quoted 
values (Detweiler,1972 ; Prec,1949). The discrepancy can be 


attributed to measurements at sites other than the aortic 
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root and/or higher flowrates than the sedentary 


6 litres/min. 


4.6 Frequency and Waveform 

The frequency of the pulsations depends on the same 
Parameters as _ the flowrate (sec.4.4). The resting rate can 
be as low as 50 beats/min. for a well-conditioned person 
(McDonald,1974) and peak at 270 beats/min. during short 
exhaustive work (Guyton,1956). An average value for an adult 
refraining from strenuous activities can be assumed at about 
70 beats/min. (Best and Taylor,1952; Altman and Ditt- 
mer,1971). The velocity waveform in the arterial system has 
been extensively analysed. As mentioned in Section 2., reli- 
able velocity waveforms can now be measured in the aortae of 
healthy people. Several publications compare graphs of 
normal and deviant velocity waveforms. For simulation, a 
waveform which appears representative of the norm has_ been 
selected from a paper by Nichols et al.(1977b) (Fig.4.5). 

The pulsation is characterized by a sharp peak during 
systole followed by a small backflow period and then minor 
OSci il omnonsmouningudiastOle sin sEne siaScenC Ingman a meac oe 
flow reversal does not occur to any significant extent 
(McDonald,1974) and the minor oscillations are largely due 
to the elastic properties of the vessel wall rather than 
backflow through the valve. Therefore, for experimental 


‘simplicity, zero flow is assumed for the complete diastole 
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period (Fig.4.6). The effects of this approximation would 


require further study. 


4,7 Dimensionless Parameters 

The agai need not be identical to the actual aortic 
arch to yield a correct flow analysis. Provided that certain 
relationships of the relevant parameters are maintained, the 
flow behaviour in the model corresponds to the physiological 
Situation. These relationships between certain dynamic and 
kinematic flow characteristics are represented by dimension- 


less parameters. 


4.7.1 Reynolds Number 

For flow through a geometrically similar model of the 
aortic arch, dynamic similarity is achieved only when the 
Reynolds number, Re, (Appendix II) is the same at every 
point in the flow at every instant in time. In other words, 
if the inlet conditions are identical in the model and 
actual aorta, setting the Re the same at two geometrically 
corresponding points ensures dynamic similarity of the flow 
everywhere in the prototype (within the previous constraints 
of inelasticity, no twist, etc). Since this experiment 
assumes uniform velocity profile at the entrance (Section 
3.), it is sufficient to assume the appropriate Re based on 
the instantaneous average velocity at the inlet of the aor- 
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In the case of steady flow, dynamic similarity is 
therefore achieved by adjusting the flowrate to yield a 
Reynolds number within physiological range. A Re of 1000 
Gornesponds Ttom thesinets positive flow in the aortic arch 
(ie.cardiac output)(Appendix II). At peak systole, Re 
attains a value circa 5000, which would be reflected in the 


model by a higher flowrate. 


4.7.2 Parameters of Periodicity 

FOC@mDULSating Wi ulOWw; met ers enOLmSUbmic Lente tTOsuSse ba eRe 
based on time-averaged flowrate since the frequency, ampli- 
tude and shape of the velocity waveform play an important 
role in the fluid behaviour. 

The unsteadiness parameter a is a non-dimensional quan- 
tity which characterizes kinematic similarities in the 
liquid motion. A typical value based on the experimental 
parameters selected earlier is a=19 (Appendix Lie 
McDonald(1974) quotes a range of 13.5<aS$16.7 at the aortic 
roou basedton a pulse of 55-72 ibeats/min.ande a “radius -sof 
1.08-1.11cem. As discussed earlier, this radius estimate is 
low. Recalculating these limits uSing R=1.45cm. yields 
Whosas 2 2.4 candles concurs with the value chosen above. 

Having specified a parameter accounting for frequency 
Similarity, it remains to define the amplitude of the wave- 
form. Decomposing the periodic flow in the cardiovascular 
system by a Fourier series, it has been shown that the first 


harmonic predominates (McDonald,1974). Therefore one form 
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previously adopted has been to approximate the fluctuations 
by a Sinusoid of a certain amplitude superimposed on a mean 
forward flow (ex. Rodkiewicz et al.,1976a). This type of 


waveform is formulated by: 


V Vee 
, = max min (4.2) 


= 


where A , the amplitude parameter, has previously been stu- 
died in the range 0< ,} <2 with the lower value representing 
steady flow and the higher value approaching an actual wave- 
Eorme (Fig74s7)' This method =does not take into account the 
Facteathat the systole sedunaticnms, onlyal/ceecot the total 
period, and that a zero flow stage is required to simulate 
diastole. However, the great advantage of this approach lies 
in its generality and ease of reproducibility. 

The experiments reported here use either steady or 
Sinusoidal flow. However, the design of a more advanced 
model (Section 5.3) adopts a characteristic velocity curve 
(Fig.4.6) and the corresponding cam profile produces the 
SAMeCMWAVELOUMewith MINOMEMOCGucactiOnsm ADDendux wll j.emennis 
method suisse comparisons with sinusoidal pulsations and 
should the differences be insignificant, the latter could be 


adopted without reservations. 
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Figure 4.1 Parts of the aorta 
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Figure 4.4 Aortic arch dimensions 
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Figure 4.5 Actual aortic pulsation(after Nichols, 1977b) 
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Figure 4.6 Design pulsation (Model B) 
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5. Apparatus and Instrumentation 
As discussed in Section 4, several simplifying assumptions 
were required to produce a feasible model. All of the tests 
carried out in this experiment make use of one model, A, 
which includes some characteristics of an actual aortic 
arch. Concurrently, a second model, 8, was designed to 
incorporate several improvements over model A and serve for 
future experiments. These models are described in detail in 


the yioltowing sections. 


See LOWRGLOCULT 

The sitlow circuit sfor, each of the two models “is /funda— 
mentally the same. A schematic and picture of the flow cir- 
cuit are shown in Figures 5.1 & 5.2. The working Bude must 
be replaced periodically as the dye used for visualization 
eventually darkens the water. To replace the fluid, fresh 
tap water is settled in an auxiliary reservoir for a period 
of several hours to prevent air bubbles from settling out in 
the model. 

A rotameter is used to roughly set the total flowrate 
inmecheCumnode! selhiesdistrpburioneolernemal Oweto thes oranches 
varies with the total flowrate through the model 
(Section 6.). To maintain a distribution within the physio- 
Tedicalerange pr raatlow Gestuictone sons séachmebrancheumustembe 
adjusted as a function of total flowrate. Orifice plates are 


used on model A while conic valves are installed on model B. 
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To monitor the flowrate in each branch, individual collec- 
tion tanks labelled T1 to T5 are used in conjunction with a 
stopwatch. These tanks in turn drain into a common reservoir 
which feeds the pump to complete the circuit. 

The. surge»tank, partially» filled) with water; acts as a 
buffer for pump fluctuations and traps air bubbles generated 
upstream. The piston apparatus is an optional attachment to 
Superimpose regular pulsations on the steady pump flow. Fin- 
oye mas e20metoot) 1ongmestratghtP sectionpamexceedingm the 
entrance length required to attain fully developed flow at 
maximum physiological Reynolds number, ensures that the flow 


conditions upstream of the model are known and consistent. 


5.2 Model A 

This model is scaled to three times actual size and its 
centreline lies in one plane (Fig.5.3). It was fabricated by 
machining two slabs of acrylic and then gluing the two sym- 
‘metrical halves together (Fig.5.4). The outlet sections are 
all positioned at the same level, about two inches above the 
model centreline. The flow distribution is controlled by a 
removable orifice plate in each outlet section. 

Thesinletesection (Fig.5.5))) ~conSistSasof an ssenlarged 
section with a transverse screen to flatten the velocity 
Duct oem rnems 1 Cnt yeclaneommode lean leiguie temcaptlU le SamsEne 
core. flow thus producing a uniform velocity protile with a 


Vetyetiim boundary layer ate che: entrancestoutnesarch. 


otes “Ue 

7 

TU ~ ov a ¢ 5 hte: ? 
5 Sr 


¥ Ae If b rl yor - i A zeq- 


. 24 
fu 4 “ve 
psy aha en e 
‘7as6) tips i AJo02 Ut 6 uf 
fi oe 
oy a 4 4 Fe Le aie 
= , 
#1 
ay fL< rh g£o.63 


i ded pot? 


Sebsciae, ay So sae fdas Pye. 


peity = 


4 
ea 
bo 
© 
sy 
J. 


en he ads a 


pon 


tev ae. Os 595. tes Se tEezT 


, 
4 i 
- @ 


chs) 


5.3 Model B 

This model was built to incorporate several improve- 
ments over model A. The principal alteration was to allow 
BOmeetnesilS° twist of the arch (Pige5. 6)tetThisswassachieved 
by splitting the model into five separate sections each of 
which can be rotated about the model centreline at the 
interfaces of sections (Fig.5.7). Hence the model can assume 
anymeangle "ore twist in@che srange sof tabouce—10muto 30>.) To 
achieve this capability however, visualization had to be 
compromised. 

Several design modifications were made to the inlet 
section of this model as well (Fig.5.8). To better assure a 
uniform velocity profile at the entrance to the arch, the 


following design factors were considered: 


ane OUT Successive screens will cause a thorough 
transfer of momentum across the inlet section to 
produce a uniform velocity profile (Appendix V). 

b. a bell-shaped entrance is provided at the arch inlet 
COMIN DLE flowsdusturbances. 

Gun Lluiders bled from) the inlets "sections sbehind = sthe 
bell-shaped entrance to streamline.the flow and pre- 


vent vortices from forming. 
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5.4 Piston 

The piston is cam-driven by a 2-Hp DC variable speed 
motor (Fig.5.9). The piston superposes fluctuations on the 
Steadvyea (Elowlesuppliced sebyeethne® pumpyaSincesthesosciitlation 
apparatus waS appropriated from a previous experimental 
set-up, the cam design required retrofitting to the existing 
apparatus. In other words, the fixed bore and limited stroke 


of the cylinder impose constraints on the cam design. 


Jem ail 

The velocity profile shown in Figure 4.6 is a simplifi- 
cation of an actual human pulsation. The backflow period and 
minor diastolic oscillations are neglected to simplify the 
cam design. It is postulated that this change would not 
Significantly affect most aspects of the flow in the model. 

Initial attempts to design a cam based on the approxi- 
mate profile yielded a pressure angle in the order of 
50°-60°. 9Since it .was found that a feasible cam profile 
could not be easily produced, a computer subroutine was 
written to produce the appropriate cam outline from input 
data containing the desired displacement curve, base circle 
radius, amount of offset and follower radius (Appendix III). 
This subprogram would calculate the maximum pressure angle 
EOuBaRGUVensprotile sandepl ctebheenesulbing contour: 

A sinusoidal pulsation was applied to model A to exa- 


mine its effects on the separation zones. The cam was 
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designed using ) =1.9 and an average Re of 1000. 

For eventual use on model 8B, a cam was designed which 
would hodtce the waveform shown in Figure 4.6. This could 
lead to studies on the effect of waveform on various flow 
phenomena in the arch. The computer program would facilitate 


the production of various other waveforms as well. 


5.6 Pulsatile Flow Monitoring 

An Annubar flowmeter was installed in the pulsatile 
flow circuit to monitor the fluctuating flow and ensure that 
the imposed oscillations are producing the desired result. 
The Annubar, a device using a calibrated Pitot tube arrange- 
ment which yields an output proportional to the flowrate, is 
connected to a DP45 Validyne pressure transducer. This 
Signal is picked up by carrier-demodulator and transmitted 
to a oscilloscope. A low-pass filter suppresses unwanted 
frequencies from the signal. 

The Annubar is specifically designed for steady flow 
but was deemed sufficiently accurate for the extremely low 
frequencies of fluctuations used in this test. However, for 
a better response, particularly at the higher frequency 
required for model B, an electromagnetic flowmeter is recom- 
mended. Some of the instruments discussed in Section 2. can 
be adapted to detect flow direction, which is advantageous 


for actual waveform simulation. 
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5.7 Flow Visualization 

The flow patterns are visualized by injecting a fine 
stream of dye through 0.5mm 1.D.tubing. The dye was alterna- 
tively introduced at various positions at the entrance to 
the arch or points along the aortic walls. At the inlet, the 
dye pipe is fixed in a ring which permits angular position- 
ing "for 180° range! ®and praesence adjustment completely 
across the tonen (CEigios10).6lFiguremes 1 tmeindreates « ssome 
peripheral injection points along the inner curve wall and 
on the top of the model. The exact non-dimensionalized loca- 
tions of these wall taps are given in Appendix IV. 

It 1S very important to maintain the densities of the 
dye and working fluid as close as possible to each other. 
Since the pump causes a temperature rise of about 6°C over a 
two-hour period of continuous operation, the water density 
continuously decreases. To compensate for this SEC Cig the 
dye reservoir is submersed in a water bath which is fed from 
Phemworkingmrlusd (Fig. 5.1 2)eeelt. isuadvaluageouUsutommalncain 
a constant dye velocity to provide stable conditions during 
visualization. By using a reservoir with a large diameter 
relative to that of the dye pipe, the loss of head in time 
becomes negligible. At the entrance to the model, the dye is 
injected with about the same momentum as the inlet flow. At 
the wall taps the dye is injected with the lowest possible 
momentum which still allows adequate visualization. 

The photographs were taken with a Canon F-1 camera 


using a 55mm lens, and optional 240mm and 450mm close-up 
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lenses, either individually or superimposed. Lighting was 
provided by one or two 200 watt tungsten lamps placed under 
the model. A single sheet of mylar was also placed between 
the lamps and the model to further diffuse the light and 
provide more uniform illumination. The black and white pho- 
tographs of the flow were taken with 400 ASA _ tungsten 
Tri-PanX and the colour pictures using 400 ASA Ektachrome. 
Although perhaps not intuitively obvious, the visual distor- 
EL One ioe Genta lathes mode sdues tO sthescunveds lanermwal Sead 
flat outer surface is less than if the outer surface were 


also curved (Roussel,1971). 
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Figure 5.3 Model A 
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Figure 5.4 Fabrication of model A 


Figure 5.5 Model A inlet section 
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Figure 5.7 Model B rotatable sections 
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Model B inlet section 
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6. Shear Stresses in Derivative Branches 

One of the important considerations in the study of blood 
flow involves the effects of wall shear stress. Several 
investigators have demonstrated some possible correlations 
between shear stresses and physiological phenomena such as 
hardening of the arteries (Fry,1968) and sites of incipient 
atheroma (Rodkiewicz,1975). However, there is no definitive 
proof of these relationships since the biochemistry of these 
problems is not yet completely understood. 

Although shear stresses are not a major issue in this 
thesis, certain results have led to speculation about them 
which would require further substantiation. Therefore, part 
of the ensuing discussion attempts simply to provide leading 


questions for future research in this area. 


6.1 Flow Partitioning 

The flow distribution in each branch as a percentage of 
total flowrate is a strong function of Re as shown in Figure 
6.1. This data was collected by choosing an orifice plate 
for each outlet which would produce the appropriate distri- 
bution at Re=1000 and then setting the flowrate to several 
Re values. These results concur with earlier data obtained 
imea Single branch model (Rodkiewicz et Val ,1976c) ) which 
demonstrated that the percent flowing through the branch 


decreases with increasing Re for steady or pulsatile flow. 
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In actual elastic arteries, the average differential 
partitioning of blood among organs is controlled by neural 
and hormonal stimulation hence the functional relationship 
is much more complex than in a rigid model. However, for a 
given physiological state, say at rest, the average resis- 
tance in each branch should remain more or less constant, 
with the instantaneous resistance primarily a function of 
flowrate, as in the plastic model. Consequently, if the 
flowrate in an actual arch were steady, the trend indicated 
in Figure 6.1 should apply. 

‘Given the pulsatile flow nature, one can speculate that 
the same effect would occur based on the instantaneous Re. 
Hence at peak systole, a larger percentage of the total 
instantaneous flowrate would pass through the descending 
aorta (branch 5) than at other instants during systole. A 
schematic of this effect is shown in Figure 6.2. 

Farthing and Peronneau (1979), having monitored the 
velocity profile in dogs' aortae, derived average flow ver- 
sus time curves for each branch. By comparing flowrates in 
the branches at various instants, it was hoped that those 
results would concur with the above assumption. However, on 
closer inspection it would appear that their plots are faul- 
ty. The peak flows in a single upper branch exeeds the peak 
flow at the inlet and the sum of the net forward flows in 
the branches exceeds by far the net flow through the inlet. 

Another study by Rodkiewicz et al.(1976a) considers the 


effect of the frequency parameter a on the flow distribution 
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with a constant net forward flow. Unfortunately there is a 
large data scatter and the results are not clear. Suffice it 


to say that these issues have not fully been clarified. 


6.2 Factors Affecting Shear Stress 

For steady flow in straight circular tubes, the wall 
Shear stress Sa is proportional to the velocity gradient 
dv/dr. As the flowrate in a given tube increases, ae 
increases also. 

Relating this to the aorta, as peak systole is 
approached, the instantaneous flowrate in all five branches 
is increasing aS 1s the wall shear stress. However, since 
the percentage flowing through each upper branch is decreas- 
ing, the shear stress 1s increased less than if the percen- 
tage had remained fixed. Therefore the varying percentage 
distribution tends to damp the wall shear stresses in _ the 
upper branches and increase them in the descending aorta. 

To get an idea of the relative magnitudes of the wall 
shear stresses in the various branches due to the net for- 


ward flow, consider the case of Poiseuille flow through a 


circular tube where: 
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where r is the radial distance measured from the central 
axis and R is the tube radius. Table 6.1 shows the magnitude 
of the wall shear stresses due to net forward flow in the 
branches, estimated by the above equation and given as a 
percentage of the inlet wall shear stress. 

Unfortunately, the effects of the instantaneous forward 
flow interact | nNon=linearly. with sthesmoulsati feu shear 
stresses. Given a Berean Sinusoidal forcing function e 
SveeecomeOnmecice DlUldgg mma eCirCUlan = D1 De mst NemWalde snear 
Stress decreases aS a increases (Fig.6.3). The relative 
magnitudes of the velocities at one point in the cycle for 
several a's are shown in Figure 6.4. 

In section 4.8.2 it was shown that a value of a=19 is a 
representative value at the aortic arch inlet. However, 
Since a iS proportional to the pipe radius R, one can define 
a different alpha for each branch as: 


2 Phranch (Ene) 


O = A 
branch Relat inlet 


The values corresponding to “ealh branch are shown in 
Table 6.1. The lower values of a found in the upper branches 
would tend to raise the maximum wall shear stress due to 
fluctuations. 

To summarize, the wall shear stress in the upper 
branches are somewhat damped by the change in flow distribu- 
tion but it would appear that this effect is completely 


overwhelmed by two other effects: the smaller a in the 
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Table 6.1 Pulsatile Shear Stress 


(branch) 
t, (inlet) 


Mean TW 
Radius (mm) 

Gre! 1.48 
Sal 2.40 
33 3.89 
ei ats) 3.89 
Ano 1.34 
Tel tad 1.48 
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branches tends to increase 1, (compared to the arch), and 
the shear stresses due to the net forward flow are much 
higher in the smaller branches. 

The above Speculations are oversimplified having 
neglected wave reflections, secondary flows, elasticity and 
other factors, but the gross effects are probably a good 
assumption and could lead to a better understanding of the 
flow phenomena. The relative influences of these factors 
cannot be predicted and would require a further experimental 


research. 
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Figure 6.2 Variation in distributiongduringecycle 
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7. Separation Zones 

For the case of steady laminar flow, a separation zone can 
be considered as a region which cannot be reached by any 
Streamline originating at ‘the model “inlet, “unless flow 
instabilities are present. For pulsatile flow, this defini- 
EVone'noe longer’ rigorousily=app! les due tromtnershnittingsor 
Streamlines with time. However the basic concept that the 
separation zone constitutes a 'bubble' or captured vortex 
distinct from the main flow remains intact. Once the loca- 
tion of a suspected separation zone has been established, 
Gvemismminjected Mdirectly—aeinto thesearca me indication or 
reversed flow and persistance of dye presence in this region 
is accepted as further evidence of the existence of a 
separation region. 

The existence and location of Separation zones in the 
aOrenerarchuangma to@derivdtivesbr anchesunasmlcugebeenvaegcon| 
tested issue. As discussed earlier, present instrumentation 
cannot detect in vivo the localized flow reversals and vor- 
tices associated with separation zones. Also, since mathe- 
matical analyses are not refined enough to predict flow 
details at this level the only recourse lies in experimental 
models. 

This report deals exclusively with two separation zones 
in the arch proper, although earlier studies have indicated 
the presence of zones in the major daughter branches as 
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on the proximal wall just distal to the entrance. The great 
majority of the tests were performed with steady flow since 
flow visualization with dye can be difficult under pulsatile 
conditions. Initial observations did not provide concrete 
evidence of reverse flow but were tantalizing nonetheless. 
In suspected separation locales, the flow would occasionally 
reverse only for an instant, leaving doubt whether the cause 
was some aioe eter rathewm than) aateature eotirptne miftow. w lit 
was Surmised, correctly, that Re=1000 was an approximate 
lower bound at which separation was initiated. As_ the 
Reynolds number was increased to 1400, the separation zones 
became more pronounced and regular. Raising the Re in the 
model to study separation can be justified on the basis that 
Re peaks near 5000 during pulsatile flow. Also the oscilla- 
tions alter the mechanism behind separation such that flow 
reversal in these zones persists even throughout diastole 
(Re=0) so that the limiting value of Re at which separation 


occurs is primarily of academic interest. 


7.1 Inner Wall Separation 

The first separation zone occurs along the inner curve 
of the arch, near the plane of symmetry and approximately 
opposite the left carotid artery (Fig.7.2). The presence of 
this region is clearly indicated in Figures 7.3(a to £) by 
Sequential photographs of flow reversal near the point of 


maximum curvature. The main flow proceeds from left to right 
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in the pictures and the dye is injected through wall tap 
number 37 (Appendix IV):* The “photographs are! taken at 5 
second intervals which, combined with the scale drawn on 
Figure 7.3a, indicates the slow progression of the reversed 
fluid. The velocity of the reversed flow is typically about 
4% of the axial core velocity. This low velocity causes some 
dye to accumulate near the injection point and as the con- 
centration increases the dye tends to diffuse slightly in 
every direction. This high concentration zone appears of the 
photographs as a small dark area around the injection point. 

The dye stream in Figure 7.3 can be easily misinter- 
preted without the benefit of a side view. However, due to 
the constrained area at the inner curve of the model, the 
camera could not be adequately positioned for a feasible 
picture of the vertical plane. Therefore, a description and 
sketches are used to complement the photographs to provide a 
complete flow picture. 

The apparent trajectory of the dye stream away from the 
innerawal lees partly ilbusonye si nmetflectmmduringm  thessbacks 
flow period ane dye is rising and following the arch wall 
quite closely. The thickness of the zone is actually about 
3-4mm. 

When the dye reaches the apex of the separation region 
most of it is entrained by the main flow and becomes 
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Depending on which streamline the dye enters, some dye may 
follow the separation bubble to its point of reattachment. 

Note? that) no recirculation is evident in the  photo- 
graphs. The dye is injected into the fluid layer nearest the 
wall which, when reaching the separation line, forms’ the 
interface between the main flow and the separation bubble. 
DUcMEO Minorsanstabllatiesmingthe, flow withismlayern econtinu- 
ously shifts slightly with the tendency to diffuse and mix 
the dye with the main flow. Injecting the dye with a higher, 
momentum to reach the inner layers of the separation bubble 
is not practicable since the dye stream completely disrupts 
the flow in the area and observations become invalid. 

Mt Seivenys Gitniculteite, ascertainmthemlocationy off the 
reattachment point since the clearance between the dye and 
the wall cannot be properly visualized. However the results 
of two other tests provide a good indication of the reat- 
tachment pointeelocationsmDyemintroduced at thesmnextes ingec— 
tion point downstream (tap #43) does not exhibit reversed 
flow so the separation does not extend that far. Dye 
injected upstream skirts the separation region and deflects 
back towards the inner wall about 0.5 centimeters past tap 
Hage (Figey.4)0.. So, although the) separationmzone, ismotedis> 
tinctly deleanated, an approximate range for this Re (1400) 
iomeErOMmmn ls oCMm.) BUDStLeam mots Lhe min jJecta onmpoln caer OmOstocm, 
downstream. In dimensionless coordinates this range extends 
from S/D=3.40 to 3.64 (Appendix IV). A schematic of the 
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widening of the zone as flow reversal progresses is probably 
accompanied by a thinning of the region although this could 
not be detected in the tests. 

Normally, flow separation is caused by an adverse pres- 
Sure gradient. If the low momentum fluid in the boundary 
layer encounters a pressure increase that it cannot over- 
come, it separates from the wall. The continuous convergence 
Oly ehelmacontic Sanchtwould tends tomsuppressmsepanata ont acma 
decreasing cross-sectional area causes the velocity to 
increase and the pressure to lower. The mechanism behind the 
flow separation at the inner wall is therefore not immedi- 
ately evident. One hypothesis which conforms with earlier 
observations attributes the inner curve separation to. the 
presence of branches. Rodkiewicz and Kalita (1978) showed 
that for a constant diameter curved tube with no_ branches 
and 400<Re<1800, separation occurs only if the Reynolds num- 
ber is large enough and the radius of curvature small 
enough. They also considered curved tubes with a side branch 
for both the cases of constant diameter and convergence of 
the main branch. In both cases separation on the inner curve 
was observed but it was not indicated at what Reynolds num- 
ber. They state however that, as expected, convergence 
decreases the size of the separation region and prevents its 
appearance for small Re. Scarton et al1.(1976) also tested a 
model geometrically similar to the arch but without side 
branches. He did not find any separation on the inner wall 
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Re=5700. Rodkiewicz(1975) also found a separation along the 
inner wall of an open channel model of the Glee uUnnere pulsas 
tile conditions. It was suggested by other researchers that 
the results were not definitive due to the difference of 
secondary flows between a channel and a closed tube. 

The @resultsPaofisMthe present experimentsin conjunction 
with the earlier works provides substantial evidence that 
the inner wall separation is caused by the existence of the 
Wan gesbronchesgori ganatinggatmr chest Oper Oremithe march... slhe 
separation is not due to the curvature of the arch alone nor 
to the influence of secondary flows, and the Reever cence Or 
the aorta is obviously not enough to prevent it. 

The presence of the branches increases the effective 
area of the arch thus increasing the pressure and giving 
rise to an adverse presSure gradient sufficiently large to 
cause separation of the slow moving fluid on the inner wall. 
The average velocity near the inner curve is lower than that 
on the side walls due to the influx of low momentum fluid by 
the secondary flow vortices. This explains why the side wall 


fluid does not separate. 


7.2 Separation on the Outer Wall 

The second separation zone occurs on the outer curve 
MOstLCaMmmOlLmmtneu li CS tempt anChem (haere rami lomsSepa hd ulon 
region is also initiated at a Reynolds number of about 1000. 


Its existence was suspected and searched for as a result of 
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earlier work by Rodkiewicz(1975) which demonstrated that 
this was also a site for atherogenesis, albeit a minor. one. 
This separation region was also evident in the open channel 
flow model of the arch. 

Figures 7.6a,b show a side view of the backflow indi- 
cated by dye alternatively injected through taps #1 and #3. 
Since the taps lie on the plane of symmetry of the model, 
the flow should divide evenly ae IngaFaquiuems/.64a mmeHowever, 
small flow inconsistencies sometimes forced the dye to fol- 
low either the upper or lower vortex (Fig.7.7) 

There is some speculation that this area is in fact not 
a separation zone but simply the site of a complex backflow 
pattern. The arguments in favor of a Separation zone are 
that no external streamline was found to enter the region; 
the location corresponds to the separation found in the open 
channel model; there is no other plausible reason for back- 
flow to occur in his area. The reasons for doubting separa- 
tion are that no vortices are evident in the region (other 
than the large secondary flow vortices); the dye does not 
reattach but is swept away (perhaps due to the fact that it 
is in the outer streamline); there is no obvious mechanism 
to induce separation at this location. With the data avail- 
abl Camiutel SemnoOtepossiblestosgainefurtherminsighteincostnas 


phenomenon pending further investigations. 
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7.3 Separation in Periodic Flow 

The observations made for the steady flow case do not 
explicitly confirm the existence of the Sepaaien zones 
under pulsatile conditions. Preliminary pulsatile tests were 
performed with a minor backflow period as in the physiolog- 
ical situation. The flow near the walls tends to reverse 
more easily than the core flow under pulsatile conditions 
Since the momentum in the boundary layer is lower. This 
eprectiers=shown sinifigures® 76.0 So although =backilow “was 
Still present in the previously defined separation regions, 
it also occurred at almost every other wall tap examined. To 
compensate for this, the input waveform was somewhat modi- 
fied so that the average flowrate remained positive through- | 
out the cycle. The amplitude varied from a minimum of Re=50 
to a maximum of Re=4000. The flow in the separation areas 
was then compared to that at several other wall positions. 

Due to the many components in the flow circuit, it was 
deemed wise to monitor the pulsations to ensure that the 
proper waveform existed at the model inlet. The output from 
the Annubar flowmeter recorded on an oscilloscope is pre- 
Sentedmine rr gqunes 7. 9a jDasthemsfirst.s picture: mindicatecsusthe 
fundamental sinusoidal input from the cam using a low-pass 
filter. The second figure represents the actual unfiltered 
waveform. The origin of the high frequency noise is not 
known but it is not equal to the blade passing frequency of 
the pump. Two plausible explanations are the noise from the 


Annubar due ‘to the pulsatile input or wave reflections from 
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the model which is just downstream from the flowmeter. 

It was found that the flow everywhere in the arch 
except in the separation regions either continued forward or 
virtually stopped as the inlet flow reached a minimum. In 
the separation zones the converse occurred. As the inlet 
flow attained its maximum the flow reversal at the wall 
became almost non-existent and as the flowrate decreased the 
fluid in the separation zone surged upstream forming a large 
backflow area. This effect was also observed in the open 
channel model by Rodkiewicz(1975). One possible explanation 
of this phenomenon is the variation in flow distribution 
with Re (Section 6.1). As the inlet Re is decreasing, a 
larger percentage of the flow is diverted through the 
branches thus effectively lowering the pressure in the arch 
in that region and causing the low momentum fluid at _ the 
inner wall to surge backwards. 

The effect of these cyclic variations is that the blood 
exerts a fluctuating stress on the arch walls, particularly 
in the separation zones. The elastic walls deform 
accordingly but eventually may lose their resiliency due to 
the large alternating shear stresses. This phenomenon may 
contribute to the development of arteriosclerosis and subse- 
quent atherogenesis as the wall metabolism is disturbed and 


fatty lesions begin to develop. 
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Figure 7.3 Progression of inner wall 
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Figure 7.3 (continued) 
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Figure 7.5 Three views of inner wall separation 
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Figure 7.8 Velocity profiles at several times(after 


McDonald, 1974) 


Figure 7.9a Filtered sinusoidal input 


Figure 7.9b Actual unfiltered signal 
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8. Flow Patterns 

The flow characteristics of the aortic arch include some 
well-studied patterns such as the large scale vortices asso- 
Clated with secondary flows and several finer points which 
even detailed observations yield grudgingly. Most of the 
testing was performed under steady flow conditions but the 
effects of pulsations were considered when possible. 

The flow patterns are examined from two points of view. 
The approximate paths of the fluid particles from the inlet 
are determined and then more detailed aspects of the flow 
are studied by injecting dye through the wall taps. The 
major characteristics of the flow are the potential core, 
the boundary layer vortices and the branch flows. 

Since the branches affect the flow a great deal, a map- 
Ping of the source of the fluid destined to pass through an 
upper branch is useful in determining their effect on the 
main arch flow. A record of the positions of the inlet dye 
pipe corresponding to a branch streamline produced the map- 
Dingeshown SineFiguress we Thewuhree! areasmilabel leds lA asmand 
4 represent the fluid which flows respectively through each 
of the three upper branches. The flows through the right 
subclavian and right carotid are conveniently lumped into 
one region (ie.1A, brachiocephalic artery) for two reasons. 
The dye has often become quite diffuse once it has attained 
tHemepiucurcacwon in the wbrachiocephalic#auteny,mpanticular hy 
Bentstaine the = boundany eiayereamakingmel Desc laicultT mato 


ascertain its ultimate destination. Also, the bifurcation 
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occurs quite far from the arch where the actual artery has 
already twisted quite a bit (Fig.4.1), hence the geometric 
Similarity of the model is poor in this location. Knowledge 
of the destination of various inlet streamlines could also 
Drove useful) for ingectingis substances) intos “the? aorta to 
enter Mas iparticulares branch. Branch?) flows. SNOWwneinmr LouUre 
Ome 

Although the boundaries between the regions are dis- 
tinct in Figure 8.1, this is not representative of the 
experimental results. When the dye is injected in the inte- 
Buona of any of the regions the results are consistent. 
However, dye injected near the boundaries of these areas 
often fluctuates between the two adjacent branches and occa- 
Sionally three. In theory, the distinction should be clear 
and the boundaries separating the areas would correspond to 
Stagnation lines between the branches (Fig.8.3). However 
minor fluctuations in the flow cause these limiting stream- 
lines to vacillate. A quick check on the data is performed 
by noting that the branch area in Figure 8.1 represents 
about 22% of the total area which corresponds to the distri- 
bution set for the branches. 

AS predicted, the boundary layer has not completely 
developed by the outlet of the model and therefore a core of 
fluids» remains © inviscid) during its passagethrough the arch 
(Fig.8.4). The slowing of the fluid in the boundary layer is 
accompanied by a velocity component inward towards the axis 


and under this influence the core fluid accelerates. By 
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injecting two dye streams side by side, it is seen that in 
this potential flow region the velocity is lower the farther 
Onesgets sfrom the-center som curvature althougn no msquantita- 
tive data was collected in this respect. The potential core 
region shrinks as the boundary layer thickens and occupies 
about 20% of the outlet area compared to 100% at the inlet. 

The area in Figure 8.1 labelled 'core' represents the 
fluid which remains in the potential core during its passage 
through the arch. Note that the whole area 1S potential at 
the inlet and that a much larger area than that marked 
NCore *sremains inviscid forvatvieast part of 1ts course, 

The remaining unlabelled region in the figure reenes 
sents the fluid which becomes part of the boundary layer or 
becomes part’ {ote@fa Ssecondarny = "flow @iwvontex ine slarges gap 
between area 1A and the outer wall indicates the amount of 
fluid which is displaced peripherally towards the inner wall 
before the first branch is reached since none of this fluid 
enters the brachiocephalic artery. 

As the boundary layer begins to grow and the fluid 
-passes into the bend, the centrifugal force is higher on the 
faster core flow which is forced outwards. Consequently the 
outer boundary layer thins as fluid is forced peripherally 
forchesioner walle Thiseetfect 1's tclearlyeeshown eine shrgure 
GeomeDymedve sin jected sintomstie boundanyslaycrgragune G.6 
shows streamlines which remain in the potential core for 
Danteao Pee Ghe ur eet ra eC LOryemondmrd ropsmeHenmentcained soymeune 


Secondary flow vortices. When observing along the arch 
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centerline facing downstream the upper vortex rotates clock- 
wise and the lower vortex counter-clockwise. This behaviour 
iS Similar to fully developed flow through a curved tube 
except, Pthat sthe, vorticess «dornotereturnetomchnerouter wall 
along the centerline. Rather they cannot penetrate the rela- 
tively high momentum potential core fluid so they remain as 
two trapped vortices along the inner wall (Fig.8.7). 

The fluid which begins along the inner wall remains 
along the inner wall as the boundary layer thickens and when 
it 1s joined at the inner curve by the secondary flow from 
the outer wall, most of it becomes entrained in the trapped 
inner wall vortices. It is quite easy to visualize the outer 
wall vortices but when the inner wall vortices first form 
they are quite small scale and difficult to observe. Initial 
observations in this area appeared only as lateral fluctua- 
tions in the dye stream, whether viewed from above or from 
the side, which was originally attributed to irregularities 
ineethes aniete flow. es ,HOWGVereme! laterals fluctuatsoOnsy Which 
appear in both the horizontal and vertical planes can be 
caused either by fluctuations in an oblique plane or by a 
small scale vortex. 

A mirror -waS mounted to attempt to photograph’ both 
views at the instant and resolve the issue. The results of 
this method were inconclusive. A second dye pipe at the 
inlet was used to inject another dye stream very close to 
the first one. These streams remained parallel for a certain 


distance and then began twisting around each other 
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indicating the presence of an inner wall vortex (Fig.8.8). 
The appearance of these vortices first occurs approximately 
opposite to the left carotid artery which is about the loca- 
tion of the inner wall separation zone. The inner wall vor- 
tices continue to grow as they proceed downstream and fill 
about half of the tube at the outlet. 

The branches have some direct effect on the potential 
core by causing a more pronounced outwards displacement than 
that which would occur in the absence of branches. This in 
turn increases the magnitude of the outer wall vortices to 
replace the outward moving inner wall fluid. 

Some observations were made under pulsatile flow condi- 
tions although the dye diffuses more easily. The vortices 
EendinemintacCte =LhiLougnout@ ase DeLlodeoteosciliation —putetne 
pitch lengthens and shortens accordingly, much as a coil 
spring would. Unfortunately, it 1s not possible to determine 
from the observations how the strength of the vortices is 
affected by the frequency or amplitude of the oscillations. 

The outer wall vortices downstream from the branches 
are less pronounced than those upstream of the branches. 
Theremaremstwo effectsewhichecontribute to@thic. Thesarcheess 
Staute ieee: Oo; COMeBOUtBOL PU SGECULVE dibten the last branch and 
as the radius of curvature increases, the centrifugal 
effects lessen. The boundary layer is also quite thin on the 
CULCTMWa UUenavVingeStaotred OVC Temabtch MUnCem ol UemsolboLavidn 
Stagnation zone. This means that the high momentum fluid is 


already near the outer wall and a significant secondary flow 
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82 
is only induced once the boundary layer becomes’ thick 
enough. 

Summarizing, the inner wall flow consists of a thick 
boundary layer and trapped inner wall vortices which begin 
Neamt news locati ones Of@amaxnmumes curvature rmeaAnee inner wale 
separation zone is also located opposite the upper branches. 
The inner wall is therefore subjected to a relatively low 
shear stress and also to localized opposing stresses at the 
separation line. The inner wall is a prime candidate for 
atherogenesis, which can be attributed to depositions 
enhanced by the slow moving blood flow as well as wall 
damage incurred by separation zone stresses. Conversely, the 
outer wall upstream of the first branch is also subject to 
atherogenesis but there the flow is characterized by a local 
ff oweneversals abdga thinnerm@boundaty lay ereshence eau nigher 
shear stress. 

Two hypotheses could explain the above observations. 
Atherogenesis could be attributed either to shear streses 
that earemtoomnighe(Fry ,1968)tor, toovlowea(Carcs ety alba 19749 
as suggested by Scarton et al (1976). These conditions could 
be aggravated by the presence of flow reversal which could 
hamper the walls' defences. The second hypothesis assumes 
that atherogenesis sites coincide with the flow reversal 
regions and do not depend directly on the main flow shear 
Stresses. This hypothesis  isispartly “substantiated | byes the 
coincidence of atherogenesis sites and separation zones in 


the branches demonstrated by Rodkiewicz(1975). 
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There was no evidence of turbulence at any time, either 
for steady or pulsatile conditions. The transition Reynolds 
number in a curved fee isa normally= higherspthan® in” fa 
Suraughe pipe yagin mthemrorden lofesRe=50008tom10,000% Thus 
depends a great deal on the curvature ratio of the tube. 
Nerem et al(1972) performing in vivo studies on dogs, made a 
rough estimate of the critical Reynolds number as a function 
of the frequency parameter a by inducing higher frequency 
pulsations wathe drugs wwerhey Sioundmithat ke 150a0. ime the 
ascending aorta and Re*250a in the descending aorta. The 
value for the ascending aorta was measured very close to the 
aortic valve where disturbances from the heart may still be 
present. In the arch itself, which is the main area of con- 
cern in this study, Re =200a can be used as a reasonable 
estimate. For a=20, the critical Reynolds number is in the 
order of 4000, so although no disturbances were observed in 
this test, the fluid may have been on the verge of transi- 
tion during the pulsatile tests. 

There is another important consideration in the _ study 
Onn mhOWeeDatterns sein nether aontlcuarch esineceGira tn ecisncums 
stances when the carotid arteries become almost completely 
occluded and neither cleaning nor medications can clear the 
lesions and restore the flow, a bypass may be installed to 
circumvent the blockage. The bypass, normally a segment of 
Vel aeiSumiSuailelny attached macomssl Cima OlUc mo UMmecem CONVEN Jeng 
heasbyes Jocation. » Often) thembypass/itseliebecomes soccliuded, 


sometimes with dire consequences. This could be caused by 
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the biochemical state of the individual's blood system. From 
the fluid dynamicists point of view however, there could 
conceivably be some advantages in selecting an optimum size, 
location and angle at which to install the bypass. The size 
of the vein used (if available) would preferably be suffi- 
cient to supply an adequate amount of blood to the deficient 
branch. The angle of attachment between the bypass and the 
aorta can affect the flow patterns and induce unwanted 
separation zones as shown by Rodkiewicz'(1973) studies of 
bifurcations. Optimizing the location of a bypass could also 
favourably affect the percentage of blood flowing through it 
as well as minimize flow disturbances. 

AS an initial guess, it was thought that locations of 
maximum pressure in the arch would be conducive to smoother 
flow into a bypass. Attempts to monitor the pressure dif- 
‘ferences between various wall points in the same plane 
resulted in the discovery that the differences were too 
small to measure with the equipment available. The taps were 
then bled in turn and the outflow collected over a long per- 
iod of time to determine the maximum pressure locations. The 
differences in hydrostatic pressure from one point to 
another were compensated for by adjusting the heights of the 
bleed tube outlets. Although there were differences in 
volume collected after long time periods, it turns out that 
the differences were not Significant and smaller than the 
errors induced by the experimental method. This test was 


therefore inconclusive. 
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To test the effectiveness of bypasses of various sizes 
and at various locations, a research program is presently 
underway to occlude some branches in the model and install 
bypasses. The study of flow patterns in and around the 
bypasses should indicate whether these parameters critically 
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Figure 8.4 Streamline in potential core 
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Figure 8.5 Dye injected into secondary flow vortex 
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Figure 8.7 Trapped inner wall vortices 
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9. Conclusions 

The major objective of this experiment, to determine whether 
Or not. a Separation zone exists at the location of maximum 
Curvature of the aortic arch, has been resolved. The separa- 
tion zone can definitely be detected in steady flow when the 
Reynolds number exceeds 1000 and in pulsatile flow it per- 
Sists throughout the cycle. The effect of the separation 
zone on disease development can be attributed to any one of 
the following: compressive stresses on the wall towards the 
Separation line; fluctuating wall stress due to continuous 
MOVeENen te moOfemthe a Separarion ine anderectuculatl on ewit han 
the Separation bubble inhibiting adequate scouring of the 
walls and transport of harmful elements away from the area. 

Earlier results (Rodkiewicz,1975) demonstrating the 
correlation between separation locations and atherosclerosis 
predilection sites in the branches and along the outer curve 
wall, left the inner curve separation as Fie main point of 
contention. Although this separation was present in the open 
channel model (Rodkiewicz, 1975) and was seen in a 
three-dimensional model,' no concrete evidence in the form 
Stee phorodtapns Was Bobtainable  siheretore apdebace@ensued? 
Since Scarton et al.(1976) demonstrated that no separation 
occurred in a curved converging tube. His claim that the 
branches do not significantly affect the flow in the arch is 
Unpounced a. aS) Orne OuteDy. thes besubLs -OLMthiseexperiment. 
The branches are in effect the sole factor determining the 


'Dr.C.M.Rodkiewicz, Univ. of Alberta; personal communication 
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existence of the inner curve separation zone. The branches 
actwas ja duiiruser ~eloweringethe velocitvein the archeinrthat 
area, thereby increasing the pressure and causing the 
adverse pressure gradient required for flow separation. 
Without the branches, a converging tube with the same curva- 
ture ratio as the arch does not show separation for the full 
physiological range of Reynolds number. 

This study provides the final link in the hypothesis 
that each atherogenesis site in the arch corresponds to a 
separation or stagnation zone in the flow. The overwhelming 
evidence of this correlation should now provide a basis for 
further biochemical studies to determine the pathological 
mechanism of atherogenesis as it relates to stresses induced 


by the blood flow patterns. 
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Appendix I : Glossary 


ACME noma ssa —CONCltiOnmymarked by | themdeposs:ingsot small 
fatty nodules on the inner walls of the arteries, and by 


degeneration of the affected areas; also such a nodule 


amcerToselerosis:; vagtnickening andslossmoteelastilerty sore the 


walls of the arteries 
atherogenesis: the primary stage of atherosclerosis 
atherosclerosis: “a “thickening ofeand loss ofselasticity in 
the inner walls of arteries, accompanied by the forma- 


tion, Of “abhecomas 


apr iMuMe war chamber lor sone mear eawhlenerece 1 Vese DloOCdm:roOnmeche 


veins and pumps it into the corresponding ventricle 
average: a weighing with respect to spatial coordinates 
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coronary artery: artery feeding blood to the heart muscle 


dvastolessthe stage: otudmletion of themventracles 


Nematocrme:. "awe thatic; Olmvolume cof packedmred’ plocdicellisato 


volume of whole blood 
minpecardivaccwexisStingeor Occurring swith un eenem tear t 
intravascular: situated or ee within a blood vessel 
MeL Ocea ne tne hiuviMow bodys Olmanmani Mal a Grne LUctmom mumans.) 
PUMmems ene cavity Of aectubUlar organ 


mean: a weighing with respect to time 


DiltccmaveUiemm=tU1 Cd) Dab temo  OOCsspascmmC Sim Uds NC CummnnOU 


suspended material 


proximal next to or néearests the point or jattachment vcr oni] 
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sclerosis: an abnormal hardening of the walls of arteries 


Sinus lof Vacsalva: “aedilation of thes aorta et thewlevel ot 


the coronary arteries 
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Subclavians thesproximaispatou, Of the Magiwance: yy ctuthe: sarm 


or forelimb 


systole: the contractiongot the heart iby whichuthegolood fis 


forced onward 


venerucie: a chamber of the heart which receives blood from 


a corresponding atrium and from which blood is forced 


into the arteries 
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Appendix II : Non-dimensional parameters 


The purpose of this section is to derive the parameters 
describing flow in a tube under pulsatile conditions. The 
constraints are that the tube is straight, rigid and circu- 


lar. Under these conditions, the flow governing equation is 


given by: 
ON an Opa y (ay a ray) Grea) 
me 1 oe aaa ole Or 


where z is the axial coordinate. The coordinates, velocity 
components and pressure term are non-dimensionalized by use 


of the following transformations: 


+ 
r=Rr Ze= oe 
at 
v=w t=it ener.) 
er 
p = (pv )p 


where the superscript plus denotes non-dimensional quanti- 
ties and L is some arbitrary axial length. 
Introducing these quantities into equation II.1 pro- 


duces the following expression after rearranging: 
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Pho. TG eHeiel Jace ay Le = = 
b. Reynolds number: Re = WO WR (11.4) 
V 


C2 HS requency soarameten:= o = Re 


SUDSUTEMCIngwiln 4mincomirepe yields: 


2 2 
2a. OV z ] op Z av 1 av " 
RS) Ae = TIES Av ¥ Re ve ie cel (11.5) 


LCS customary to adopt the definition of a given in equa- 
tion IJI1.4 (even though it produces the parameter squared in 
II.5) since a shows up again in the solution to equation 
ely ae 

Equation II.5 can be applied to the case of pulsatile flow 
with zero net flowrate where the amplitude of the velocity 
wave is governed by the amplitude and frequency of the forc- 
has GieebaVeneallerg op . In a more realistic application for blood 
Elowestuaqvess the forcingmruncti on | cONnSisStS mOL aseconstant 
pressure gradient term superimposed by a fluctuating com- 
ponent. This produces a net flowrate component and an oscil- 
Vating flow component... Thesrelative amplitudes otsthesescom- 


ponents can be represented by a non-dimensional amplitude 


parameter (Kuchar and Scala,1968): 
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where \X represents the ratio of peak to peak velocity to 
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the steady velocity component. Of course if Vain =2 chen 


becomes: 


‘2 u CL eee) 


Heretofore, ,} hasS primarily been used in applying a 
Sinusoidal pulsation but is equally valid in scaling an 
actual pulsation to model proportions. In_ other words, to 
Simulate an actual pulsation with a cam and piston arrange- 
ment requires input of the waveform to the design program 
DUDMESCa ling» «Cane Stillmebe: performed. withetne sald otatne 
amplitude parameter. 

The dimensionless length parameter is satisfied by 
reproducing the actual geometric eens in the model. The 
Reynolds number is based on the mean inlet velocity which of 
course depends on the cardiac output so a representative 
value has been selected. Finally the frequency parameter 
a=19 as an experimental choice is based on a heart rate of 


70 beats/minute. 
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Appendix III ;: Cam Design 


It is difficult to produce a feasible cam which would ade- 
quately simulate the heart pulsations. Therefore, the cams 
used in this experiment were designed through a versatile 
computer program. The program actually consists of two dis- 
tinct subprograms. One subroutine was specifically designed 
for this experiment to generate a displacement curve from a 
representative velocity waveform of flow in the aorta. This 
displacement curve is then inputted to the cam design  sub- 
routine. The second subroutine produces the profile of a 
disk cam with reciprocating follower from input data of the 
corresponding displacement curve. This routine 1S very gen- 
eral and accepts a wide range of input parameters. 

The input to the displacement subroutine is a series of 
datappoints flrompaneaveragemvelocilyavsmcimestrace OL ler low 
= an actual aorta. To reproduce a similar waveform in the 
model, the scaling of the velocity curve must be modified so 
that the non-dimensional parameters a and Re are both satis-~- 
fied. 

The vertical scale is determined by the Reynolds num- 


ber: 
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where subscripts 1 and 2 refer to model and aorta respec- 
tively. ~Since (the. cam ssuperimposes a fluctuation over the 
steady flow from the pump, this must be accounted for by 


subtracting the steady velocity component from the instan- 


taneous velocity: 


] 
ery aa: J v,(t)dt (Ol Ee 


Fonal lye Chepeverticaleescale OL ms tnemeve Oct Ver curveymiS 
adjusted to allow for the difference between the piston bore 
and the diameter of the model inlet. 

The horizontal or time scale must be determined from 
the period of the Puc tuatwon which is inversely 
proportional to the frequency. The frequency of the flow in 


the model must be set to satisfy a: 


9 (ell 1e3}) 


Once the scaling is completed, the subprogram inte- 
Gratesmaticmivelocity sscUnvVcisbyms liNpSoOnmecerulemtomy: eldagthe 
Corresponding displacement protile sl tecanebesnoveds that tie 
vertical scale is proportional to the viscosity ratio and 
the horizontal scale is inversely proportional to this 
ratio. Therefore, when integrating, the viscosity ratio term 


drops out of the equation. In other words, the cam profile 
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does not depend in any way on the working fluid being used 
in the model. The displacement curve corresponding to the 
VelLOCirye protile from Figure 4.6 15 SnOwleineh laure itil 
The cam subroutine then reads in data points from the 
displacement curve, asks for a base circle radius and con- 
verts the horizontal scale. to degrees of rotation, bringing 
the design to a stage where the pitch profile can be 
plotted. Due to the sudden increase in velocity at the 
beginning of systole, the pressure angle on a centered cam 
PSemuchetoo high for they cam to be practicable, “Therefore, 
the program was expanded to allow for offsetting the cam to 
reduce the pressure angle when required. Finally, using 
Simple calculus, the program determines the inward normal to 
Ene witch protrie at each data point, marks off a distance 
along the normal equal to the follower radius (as specified 
byMenem programmer), and 1s sGhUSm Ins cm U0S it Onm iCOpm pL Ocemcile 
complete cam profile and quote the maximum pressure angle 
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UicweLOuN tle wavelOrM ineaneactUalsaolcamaremonOwim nl Pian 
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Appendix IV : Injection Points 


The locations of the dye injection points must be accurately 
specified sovthat»f£lowlcharacteristics adisplayedeat «Various 
locations can be integrated into a single complete flow pat- 
tern. Since thesmodel is three times) actuals @size™ “to is 
advantageous to use a non-dimensionalized coordinate system 
to denote the infusion sites. Two coordinates S and 6, are 
used where: 

S=sdistancesalong the centreline lof | themecortic, arch 

from the inlet plane 

® = angle (in radians) between a line OP normal to the 

centreline and passing through the injection point, and 

hvne, OH whien= 1s the outwardenormal torSeavong the honi— 

zontal plane of symmetry (Fig.IV.1) 
The downstream position is then non-dimensionalized by 
ivi tn Gm Yeecicmy nleleCiame le) Bote nie wc GCil neta crea o 
indicates the numbering scheme for the wall taps and Table 
IV.1 presents the coordinates of =each wall imjection point. 
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Table vive 1. Dye Tap Locations 
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Figure IV.2 Dye pipe numbering 
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Appendix V : Inlet Section Design 


Asm the velocity profile at the entrance tomtne aortic arch 
plays such a dominant role in the development of downstream 
flow patterns, it is imperative that the Dae section pro- 
vide a means for complete momentum transfer across the chan- 
nel. This is achieved by a_ series of polyester screens 
spaced evenly along the entrance pipe. The polyester 
material is precision woven and inert in water, therefore 
ideal for this application. The relevant criteria used to 
determine the appropriate mesh size and screen placement 

are: 
a. the porosity of the screen should be about 50% anda 
range of 40-60% is acceptable, where porosity 1s 


defined as: 
3 d 
ot - iy (V.1) 


where d and W are the wire diameter and spacing 
respectively; 

b. the screens should be placed ae least 200 wire dia- 
Metenseapalt. ) Witheaboutm O00 swWinercdlamecers@amver 
the final screen; 

c. four screens should provide adequate transfer of 
momentum. 


These rules of thumb are commonly used in experimental 
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practice* and have proven suitable in a wide range of appli- 
GatwonS ms (neGatirSstast oUlattony fOMlOovGmELometnemrac teria tems 
large porosity will not effectively transfer momentum, and a 
low porosity will cause a large pressure drop and unaccept- 
able jetting effects. The second empirical condition allows 
enough space for the small scale eddies to damp somewhat 
before the following screen and,the final spacing is recom- 
mended to completely dissipate any turbulence caused by the 
mesh. However, since it 1s also desirable to inhibit boun- 
dagyeleyere CuOwthmupStreampotmrne arch inlet, a compromise 
of 400 wire diameters after the last screen was settled 
UDOn a fainavly ss thesaddi tional mixingrachlevecebyeusing amore 
than four screens is usually not justified for the extra 
space and pressure drop required. 

To satisfy these requirements, a screen with a thread 
diameter of. 0.018 inches: Wand) Meshepcounty inch=21x21swas 
Selected se Ihis gives as pOrosl uy Of B= 56. sand anmwaxiale spaced 


of =2.4 inches. between screens. 
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